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We explore the possibility of describing the main transport properties of a granular gas by means 
of a model consisting of elastic hard spheres under the action of a drag force that mimics the inelastic 
cooling of the granular gas. Direct Monte Carlo simulations of the Boltzmann equation show a good 
agreement between the results for a gas of inelastic hard spheres and those for a gas of driven elastic 
hard spheres in the simple shear flow state. This approximate equivalence between both systems is 
exploited to extend known kinetic models for elastic collisions to the inelastic case. 



I. INTRODUCTION 



The simplest model of a granular fluid in the rapid flow regime consists of a gas of (smooth) inelastic hard spheres 
(IHS) characterized by a constant coefficient of normal restitution a The inelasticity of collisions produces a 
decrease of the mean kinetic energy (or granular temperature) with a cooling rate C oc 1 — a 2 . Interestingly, a cooling 
effect can also be generated in a gas of elastic hard spheres (EHS) by the application of an effective drag force 
with a friction coefficient |£. At a macroscopic level of description, the hydrodynamic balance equations of mass, 
momentum, and energy for the IHS gas are (formally) identical to those for the frictional EHS gas. However, the 
microscopic dynamics is physically quite different in both systems: in the IHS gas (a) the particles move freely between 
two successive collisions but (b) each colliding pair loses energy upon collision; in the EHS case (a) the particles lose 
energy between collisions due to the action of the drag force but (b) energy is conserved by collisions. This implies 
that during a certain small time step, only the small fraction of colliding particles are responsible for the cooling of 
the system in the IHS case, whereas all the particles contribute to the cooling in the EHS case. Therefore, there 
is no reason in principle to expect that the relevant physical properties (e.g., the velocity distribution function) are 
similar for IHS and frictional EHS under the same conditions. For instance, in the so-called homogeneous cooling 
state the solutions to the respective Boltzmann equations for IHS and EHS differ: while the distribution function is a 
(time-dependent) Gaussian for EHS |2|, deviations from a Gaussian (as exemplified by a nonzero kurtosis and by an 
overpopulated high energy tail) are present in the case of IHS Q. Notwithstanding this, the differences between the 
homogeneous solutions for IHS and EHS are not quantitatively important in the domain of thermal velocities and so 
it is still possible that both systems exhibit comparable departures from equilibrium in inhomogeneous states where 
transport of momentum and/or energy is the relevant phenomenon. The investigation of this possibility is the main 
aim of this work. 



II. MODEL OF DRIVEN ELASTIC HARD SPHERES 



The Boltzmann equation for a gas of inelastic hard spheres (IHS) is 0,0 

(a t +v-V)/(v) = a*- 1 J dvi I d^e(g-^)(g-^)[a- 2 /(v')/(vi)-/(v)/(v 1 )] 

= J (a) [/,/]■ (1) 

In this expression /(r, v; t) is the one-particle distribution function, a is the diameter of a sphere, d is the dimensionality 
of the system, is the Heaviside step function, <r is a unit vector directed along the centers of the two colliding spheres 
at contact, g = v — vi is the relative velocity, and a is the coefficient of normal restitution. The precollisional or 
restituting velocities v' and are given by 

, ^ + ol ^ l + a, 

v =v 2a~(S' <T ) cr ' v i = v i + -^(g ' <T ) (T - (2) 
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The collision operator for elastic hard spheres (EHS), J^[f, /], is obtained from eqs. QJ and J5J by setting a = 1. 

The first d + 2 moments of the distribution function define the number density n, the noncquilibrium flow velocity 
u and the granular temperature T. The most important properties of J( a > [/, /] are those that determine the form of 
the macroscopic balance equations for mass, momentum and energy, namely 

J dv{l,v,my 2 } [/,/] = {0,0, -dnTQ, (3) 

where m is the mass of a particle, V(r,i) = v u(r,i) is the peculiar velocity and ((r,t) is the cooling rate due to 
the inelasticity of the collisions. Although it is a nonlinear functional of the distribution function / 5, 6], a simple 
estimate £ w £o is obtained by replacing the actual distribution function / by the local equilibrium distribution 
fo = n(m/2irT) d / 2 exp{-mV 2 /2T). The result is [10 

,rs , N rf + 2 /, 2^ Stt^-D/V*- 1 (T\ 1/2 

C (r,t)=Mr,t)— r (l-a 2 ), v = , , w W M ~ • ( 4 ) 
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4d v " u (d + 2)T(d/2) \m 
According to the arguments of section ^ our model consists of the replacement 

J {a) [fJ] - /?(a)J (1) [/,/] + \Co(a)^ ■ (V/) 



where /3(a) is a positive constant to be determined. The model J5J complies with the properties ©, except that the 
true cooling rate is replaced by the local equilibrium estimate Q • According to eq. ijjjj , the gas of inelastic hard spheres 
is replaced by an "equivalent" gas of elastic hard spheres subjected to the action of a drag force F<j rag = — (m£ /2)V 
proportional to the peculiar velocity with a friction constant that depends on the local density and temperature. This 
drag force mimics the cooling effect due to dissipative collisions in the underlying granular system. The parameter (3 
accounts for the fact that, in principle, the gas of EHS that more efficiently succeeds in capturing the main properties 
of the granular gas is made of particles with a diameter a' — that does not necessarily coincide with the 

diameter a of the inelastic spheres. Alternatively, we can view (3 as a correction factor to modify the collision rate 
of the equivalent system of EHS. A comparison between the transport coefficients of IHS 0, and those of the 
"equivalent" EHS yj suggests the choice 

/3(a) = -(! + «). (6) 



III. SIMULATIONS 



In order to test the model (|5jl in inhomogeneous states far from equilibrium, we have carried out computer simula- 
tions of the Boltzmann equation by means of the DSMC method [l(| in both systems (IHS and EHS) for the simple 
shear flow problem. In the simple shear flow |n| , the gas is enclosed between two infinite parallel plates located at 
y = ±L/2 and moving with velocities ±C7/2 along the a;-axis. When a particle crosses one of the plates it is reentered 
through the opposite plate by applying the standard Lees-Edwards boundary conditions 01 ■ This produces a viscous 
heating effect that tends to increase the temperature of the system, whereas the inelastic cooling (in the IHS system) 
or the drag force (in the EHS system) tend to decrease the temperature. Eventually, a noncquilibrium steady state 
(NESS) is reached when both effects cancel each other. If the size of the system is large enough as to avoid clustering 
effects, the NESS is characterized by uniform density and temperature, and a linear velocity profile u x (y) — ay, where 
a = U/L is the constant shear rate. 

As a test case, we have considered three-dimensional systems with L — 2.5A, where A = (s/2'Kfia 2 )^ 1 is the average 
mean free path of the IHS gas (n being the average density), and U = I0v , where vq — \J 2T / m is the initial thermal 
velocity (T being the initial temperature). The shear rate is then a = Atq , where r = X/vq is the initial mean free 
time of the IHS gas. The coefficient of restitution for the IHS gas has been taken asa = 0.9. As for the EHS gas, its 
collision rate has been reduced by a factor (3 = 0.95, in agreement with eq. ©. In both cases the viscous heating effect 
dominates during the transient regime until the NESS is reached. In the simulations we have considered a number 
N = I0 4 of simulated particles, a width layer 5L = 0.05A and a time step St — 10 -3 to\/To/T. 

Figuren^a) shows the velocity and temperature profiles at times t/ro = 0.13, 0.5, 1, 1.5, 2 for both systems, starting 
from an initial condition of total equilibrium. By time t — 2tq (what corresponds to about 2.4 collisions per inelastic 
particle) the velocity profile is practically linear and the temperature profile is almost uniform. However, the global 
temperature keeps growing in time until it reaches a stationary value T ~ 146Tb for t > 10r (i.e. after about 48 
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collisions per inelastic particle). The time evolution of T/Tq and of —P xy /nTo (where P xy is the shear stress) is 
shown in fig. ^b) for an initial condition of local equilibrium, so that the system is initially prepared with a linear 
velocity profile. The NESS value of the shear stress represents an effective shear viscosity about 14% smaller than 
the Navier-Stokes value corresponding to a — 0.9. Figure ^ along with a more comprehensive comparison that 
will be reported elsewhere 0, shows that the equivalent EHS system succeeds in capturing the main nonequilibrium 
transport properties of the underlying IHS system. 
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FIG. 1: (a) Velocity and temperature profiles, (b) Time evolution of the temperature and the shear stress 



IV. KINETIC MODELING 



A. BGK and ES models 



The mapping IHS^EHS allows one to take advantage of the existence of simple kinetic models for EHS to extend 
them straightforwardly to IHS. For instance, consider the so-called ellipsoidal statistical (ES) model [Tll[T^ | 



where 



J (1) [/,/]--^o(l 

/.(v) = n (^f 2 (det R)-/ 2 exp (-^R- : w) , R = ^ (p. - eP) , 



(7) 



(8) 
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P being the pressure tensor. The parameter e < dT 1 can be freely chosen. The choice e = d _1 makes the ES model 
reproduce the correct value Pr = (d — l)/d of the Prandtl number. The simplest choice, however, is e = 0, in which 
case R = pi, the reference function f e = fa becomes the local equilibrium distribution and the ES model reduces to 
the well-known BGK model. 

In the spirit of (JSJ, the extension of the ES model to IHS is 

J (a) [f, /] - -/3(aH(l - e)(f - h) + \Co(a)-^ ■ (V/) . (9) 

In particular, setting e = we get a simplified version of the BGK-like model for IHS that was proposed in ref. 
This BGK model is easy to solve for the simple shear flow problem considered in section IIIII The thin solid lines 
in fig. ^b) represent the evolution of the temperature and the shear stress according to such a solution. As can be 
observed, the BGK solution exhibits an excellent agreement with the DSMC results for both IHS and EHS. 



B. Mixtures 



The same idea behind @ can be extended to a multi-component granular gas [9j . In the special case where all the 
species have the same flow velocity (uj = u), our model becomes 

^hfiJA^PiM^ih'fA + l^M^ ■ K v - u *)/i] ( 10 ) 

with 

where Ti is the granular temperature of species i and (iji = (1 + mi/rrij)^ 1 . The model (|1()J) preserves the first 
d + 2 collision integrals of IHS in the leading Sonine approximation with = u. The important point is that the 
approximation (|10fl allows one to transfer any given kinetic model 

jSh^JA^K^ (12) 
for elastic mixtures [TlT | into an equivalent model for inelastic mixtures: 

4T w) [/- fi\ - K t 3) - X - ± y Lk ^ + C f^ ■ [(v - m) /*] ■ (13) 



V. CONCLUSION 



In summary, we have shown that the nonequilibrium transport properties of a Boltzmann gas of inelastic hard 
spheres can be satisfactorily captured by an equivalent gas of elastic hard spheres driven by a dissipative drag force 
Fdrag = — (m£o/2)V, where Co( a ) is the (local equilibrium) cooling rate. Besides, the elastic particles must reduce 
their collision rate by a factor (3(a) w |(l + a) in order to "disguise" as a granular gas. While the "equivalent" system 
of EHS does not retain finer details of the true IHS gas (e.g., high energy tails, velocity correlations, . . . ), it is able 
to account for those phenomena (e.g., inelastic clustering) that can be described at a hydrodynamic level. Finally, we 
have exploited the possibility of reverting the mapping IHS^EHS to construct kinetic models for granular gases as 
natural extensions of known kinetic models originally proposed for elastic particles. 
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